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"Method and apparatus for magnetic field detection" 
Technical Field 

The present invention relates to a method and apparatus for 
measurement of a magnetic field, and in particular relates to rotating a 
5 magnetic sensor relative to the magnetic field to be measured. 

Background Art 

A gradiometer is a device for measuring the difference between two 
signals, with an aim to facilitate rejection of common mode noise signals and 
10 improved reduction in errors that arise from misorientation of the measuring 
instrument. 

The ability of a gradiometer to measure small magnetic field gradients in 
the presence of a large background magnetic field is determined by its ability to 
reject large common mode signals, and the intrinsic noise and drift 

15 performance of the field detectors employed. The first of these characteristics 
is called the common mode rejection ratio (CMRR). In conventional 
gradiometers the CMRR is determined by the precision with which the field 
sensors comprising the gradiometer are balanced. That is, the CMRR is 
determined by the relative sensitivities of the field sensors, and the accuracy to 

20 which the field sensors can be positioned in precisely parallel planes. This 
fundamental nexus between the precision with which conventional 
gradiometers can be manufactured and the ultimate CMRR performance has 
until now limited the CMRR to figures of the order of 10" 4 in practical 
instruments. 

25 Further, the DC and low frequency accuracy of state of the art 

gradiometers is limited by the DC and drift performance of the field detectors 
employed, and is also limited by the low frequency noise performance of the 
field detectors. Furthermore, the achievable noise performance can be 
dependent upon the magnitude of the background homogeneous field; being 

30 determined by microphonics which arise from vibrations causing randomly 
varying misalignment of the axes of symmetry of the two field sensors. For 
example, in the case of magnetic gradiometry two appropriate magnetic 
sensors are the flux gate and the superconducting quantum interference device 
(SQUID). Flux gates have considerable drift and noise in comparison to 

35 SQUIDs, but state of the art SQUID systems are incapable of providing data 
about the absolute value of the magnetic field or gradient. 
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The inability of state of the art SQUlDs to measure the absolute value of 
the magnetic field or gradient arises for the following reasons. The radio 
frequency (rf) SQUID is essentially a ring of superconducting material, the ring 
being interrupted by a Josephson Junction. When the superconducting ring is 

5 energised by an inductively coupled resonant rf-oscillator, tunnelling of 
electrons takes place at the junction and a periodic signal, being a function of 
flux through the ring, can be detected across the junction. The periodic signal 
is substantially a triangular waveform, usually having a period (AB) in the order 
of a nanotesla. Therefore, in order to yield a sensitivity in the femtotesla range, 

10 the SQUID is operated in a nulling bridge mode, or flux locked loop (FLL) 
mode. In this mode, magnetic flux is fed back to the SQUID so as to cause the 
output voltage to remain relatively constant. The feedback voltage, being 
proportional to the difference between the applied flux and the quiescent flux 
level, gives a highly accurate measurement of relative magnetic flux. The 

15 feedback voltage V can therefore be written as 

V = M (AeffB + u) (1) 

where 

M is a constant in a specific SQUID system; 
A eff is the effective area of the SQUID; 
20 B is the applied magnetic field; and 

u is the quiescent flux. 

However, the quiescent flux u is unknown. Thus, conventional rf- 
SQUIDs provide only relative measurements of magnetic field, up to an offset 
(which can usually be held constant for the duration of the measurement 

25 procedure) which is an integer multiple of a field value determined by the 
sensitivity of the SQUID and the fundamental flux quantum. The SQUID, 
therefore, is not an absolute field detector and state of the art gradiometers 
cannot provide absolute value measurements of the gradient components. 

A further problem arises in tensor gradiometers. The gradient of a field 

30 in 3 dimensions is a second rank tensor, and accordingly has 9 components. 
In the case of magnetic fields Maxwell's equations impose restrictions which 
imply that only 5 of the gradient components are independent. To determine 
the gradient tensor of a magnetic field it is therefore necessary to measure all 5 
of these components. State of the art tensor gradiometers achieve this by 

35 deploying a suite of at least 5 field sensors or planar gradiometers, or a 
combination of field sensors and planar gradiometers. 
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However, a problem which arises from the use of a plurality of field 
sensors or gradiometers to measure the gradient components, is that the 
sensors or gradiometers distort the field in their neighbourhood. Consequently, 
each sensor or gradiometer adversely affects the measurement of the field by 
5 the other sensors or gradiometers which are deployed to measure other 
gradient components. Such a "cross-talk" effect limits the accuracy of the 
gradient measurement. 

Any discussion of documents, acts, materials, devices, articles or the like 
which has been included in the present specification is solely for the purpose of 
10 providing a context for the present invention. It is not to be taken as an 
admission that any or all of these matters form part of the prior art base or were 
common general knowledge in the field relevant to the present invention as it 
existed before the priority date of each claim of this application. 

Throughout this specification the word "comprise", or variations such as 
15 "comprises" or "comprising", will be understood to imply the inclusion of a 
stated element, integer or step, or group of elements, integers or steps, but not 
the exclusion of any other element, integer or step, or group of elements, 
integers or steps. 

Throughout this specification, the terms 'superconducting material 1 , 
20 'superconducting device' and the like are used to refer to a material or device 
which, in a certain state and at a certain temperature, is capable of exhibiting 
superconductivity. The use of such terms does not imply that the material or 
device exhibits superconductivity in all states or at all temperatures. 

25 Summary of the Invention 

According to a first aspect, the present invention provides a method of 
field detection comprising: 

providing a gradiometer in a field, the gradiometer having at least first 
and second field vector sensors connected in a differencing arrangement; and 
30 controllably altering a position of the at least first and second field vector 

sensors relative to the field during operation of the gradiometer. 

According to a second aspect the present invention provides a field 
detection device comprising: 

a gradiometer, the gradiometer having at least first and second field 
35 vector sensors connected in a differencing arrangement; and 
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means for controllably altering the position of the at least first and second 
field vector sensors relative to a field during operation of the gradiometer. 

In preferred embodiments of the present invention, the at least first and 
second field vector sensors of the gradiometer are substantially axially aligned, 
5 thus forming an axial gradiometer. That is, the at least first and second field 
vector sensors are preferably positioned relative to each other such that first 
and second field vectors sensed by the respective first and second field vector 
sensors are substantially co-axial. However, the present invention may also be 
applied in respect of other devices such as planar gradiometers. 

10 In particularly preferred embodiments of the present invention, the 

position of the at least first and second field vector sensors is controllably 
altered by rotating the at least first and second field vectors about an axis of 
rotation. Such rotation may be continuous, as in the case of a motor driving the 
gradiometer at a substantially constant revolutional velocity about the axis of 

15 rotation. Alternatively, the rotation may be piecewise about the axis of rotation, 
as in the case of manual rotation of the gradiometer through a given angle 
about the axis of rotation. Measurements may be obtained from the 
gradiometer before, during, and/or after such a manual rotation. 

In preferred embodiments of the invention, the axis of rotation is 

20 . positioned substantially perpendicular to the co-axial first and second field 
vectors. The axis of rotation preferably passes substantially between the first 
and second field vector sensors, and substantially equidistant from the first and 
second field vector sensors. 

In accordance with preferred embodiments of the present invention, the 

25 gradiometer is made to rotate so that the direction of the sensitive direction of 
the field vector sensors or their sensitivities is changed in a known way. 

In accordance with preferred embodiments of the invention, the 
gradiometer may be an axial magnetic gradiometer, and may be formed from a 
pair of SQUIDs, flux gates or superconducting pick-up loops with their 

30 sensitivity vectors lying substantially in a nominal x-y plane. In preferred 
embodiments, such an axial gradiometer is rotated about a nominal z-axis to 
form the rotating gradiometer. 

While the use of an axial gradiometer may be favoured in 
implementations of the present invention because of its intrinsic noise 

35 immunity, the advantages of rotating a gradiometer are not limited to the axial 
configuration. 
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It has been realised by the present inventors that a rotating axial 
magnetic gradiometer returns information about the value of the gxy component 
and linear combinations of the values of the gxx and gyy components of a field 
gradient tensor, in addition to information about the B x and B y field components. 
5 In preferred embodiments of the invention, such information is retrieved by 
transforming the voltage waveform from the field vector sensors into the Fourier 
domain. 

Furthermore, it has been realised that it is possible to obtain the 
complete magnetic gradient tensor by deploying three such rotating 
10 gradiometers, configured in such a way that their respective axes of rotation 
are not parallel. 

Accordingly, in a third aspect, the present invention provides a method of 
obtaining a complete magnetic gradient tensor of a magnetic field, the method 
comprising: 

15 providing at least three axial gradiometers in the magnetic field such that 

an axis of each axial gradiometer is not parallel to an axis of any other one of 
the at least three gradiometers; and 

controllably altering a position of each of the at least three gradiometers 
relative to the magnetic field during operation of the gradiometer. 

20 Preferably the step of controllably altering comprises rotating each of the 

at least three axial gradiometers about a respective axis of rotation of that 
gradiometer. The axis of rotation of each gradiometer is preferably not parallel 
to an axis of rotation of any other one of the at least three gradiometers. 

Similarly, in a fourth aspect, the present invention provides a device for 

25 obtaining a complete magnetic gradient tensor of a magnetic field, the device 
comprising: 

at least three axial gradiometers positioned such that an axis of each 
axial gradiometer is not parallel to an axis of any other one of the at least three 
gradiometers; 

30 means for controllably altering a position of each of the at least three 

gradiometers relative to the magnetic field during operation of the gradiometer. 

Preferably the means for controllably altering comprises means for 
rotating each of the at least three axial gradiometers about a respective axis of 
rotation of that gradiometer. The axis of rotation of each gradiometer is 

35 preferably not parallel to an axis of rotation of any other one of the at least 
three gradiometers. 
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In preferred embodiments of the third and fourth aspects of the invention, 
the at least three gradiometers are rotated at a substantially constant velocity 
about their respective axes of rotation. Such embodiments produce data 
suitable for conversion to the Fourier domain. 

Embodiments of the third and fourth aspects of the invention effectively 
provide a single measurement of the true value of each of the three off- 
diagonal components of the gradient tensor, three measurements of the 
diagonal components and two measurements of each of the field components. 
Because more than the minimum of five tensor components are obtained, the 
use of just three such rotating gradiometers in fact provides a redundant 
measurement. This redundancy permits a least-squares best fit solution, 
further improving the quality of data obtained. Any DC offsets such as low 
frequency drift in flux gates or the fixed offset of a SQUID, for example, can 
also be determined and monitored to provide information about the operating 
conditions of the detectors. Such offsets can also be separated from gradient 
data in the Fourier domain. 

In preferred embodiments of the third and fourth aspects of the invention, 
data obtained is converted to the Fourier domain. In the Fourier domain the 
gradient information is distinguishable from the field information because the 
signals due to the gradient fields are centred at twice the rotation frequency 
while those due to the fields are centred at the rotation frequency. 

Further, information about the gxy component of the gradient tensor may 
be distinguished from that due the diagonal gradient components, even at the 
same frequency, because one involves taking the Fourier transform of a cosine 
waveform, which is real and even, while the other involves the Fourier 
transform of a sine wave, which is odd and imaginary. The Fourier 
transformation may be implemented by means of a conventional fast Fourier 
transform algorithm or even simple discrete implementation via a coarse 
sampling of say 8 data points per cycle of the rotation. 

To first-order any mismatches in the sensitivities of the field sensors or 
their relative alignment gives rise to contributions only at the fundamental 
frequency (discussed in more detail in the following) and therefore such 
mismatches do not determine the ultimate CMRR of the gradiometer. The final 
CMRR performance is determined by the resolution arid precision of the data 
handling and the ultimate noise floor of the detectors. Accordingly, preferred 
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embodiments of the present invention may provide for a CMRR improved by 
several orders of magnitude over that of state of the art instruments. 

Further, since the external gradient and field are modulated by the 
rotation, but DC offsets, drift or very low frequency noise in the detectors are 

5 not, the rotation provided by the present invention facilitates separation of the 
these effects in the Fourier domain so that the DC and low frequency 
performance of the instrument is not limited by the performance of the 
detectors in these respects. 

Furthermore, in the case of magnetic gradiometry, if the gradiometer is to 

10 be employed in a moving platform such as an aircraft, for example, a 
suspension system is generally used to isolate the instrument from vibrations 
which would otherwise induce significant noise by moving the instrument in the 
earth's magnetic field. However, in previous systems, such suspension 
systems have proven to be extremely difficult to engineer as significant 

15 isolation is required at sub 1 Hz frequencies. However, the rotating 
gradiometer provided by present invention modulates low frequency signals 
with the rotation frequency, thus facilitating the use of suspension systems 
which then only need to provide effective vibration isolation near the 
fundamental and first overtone of the rotation frequency. Such suspension 

20 systems are typically significantly easier to engineer than sub-1Hz suspension 
systems. 

In preferred embodiments of the third and fourth aspects of the invention, 
the at least three gradiometers are rotated at differing frequencies. Such 
embodiments of the invention facilitate separation of the frequency of operation 
25 of each gradiometer, and further facilitate a reduction in "cross-talk" between 
the respective gradiometers, of the type set out in the preceding discussion. As 
such cross-talk effects will be occurring at differing frequencies, the effect of 
one gradiometer operation on another can be filtered out or removed in the 
frequency domain. 

30 Accordingly, preferred embodiments of the present invention may yield 

significant advantages over previous field detection devices and methods in a 
number of ways. First, rotation of the gradiometer breaks the nexus between 
engineering tolerances (in positioning two sensors in precisely parallel planes) 
and the ultimate CMRR performance. Further, providing three rotating 

35 gradiometers in accordance with embodiments of the third and fourth aspects 
of the present invention yields measurements of all five components of the 
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gradient tensor and all three the components of the field, with redundancy in 
the data, all by using only 3 field detectors. Additionally, the rotation of the 
gradiometer breaks the dependence of the DC performance of the instrument 
on the DC performance of the detectors. Embodiments of the present invention 
also facilitate absolute value measurements of both the gradient and field 
components using SQUIDs. Embodiments of the present invention may also 
provide real time information of the condition of the SQUID in operation, for 
example as to whether or not the SQUID has remained in lock. As noted in the 
preceding discussions, previous SQUID devices have only been capable of 
making relative field measurements, not absolute field measurements. 
Furthermore, embodiments of the present invention may improve the sensitivity 
of field detection by shifting the working point of the detectors to a higher 
frequency, where the noise performance of the detectors and the vibration 
isolation afforded by suspension systems is superior. Yet further, certain 
embodiments of the present invention can serve to reduce the effects of cross- 
talk by use of distinguishable rotation frequencies. 

Brief Description of the Drawings 

By way of example only, preferred embodiments of the invention will be 
described with reference to the accompanying drawings, in which: 

Figure 1a illustrates a tow temperature superconducting gradiometer in 
accordance with the present invention; 

Figure 1b illustrates a high temperature superconducting gradiometer in 
accordance with the present invention; and 

Figure 2 illustrates a three-gradiometer magnetic field detection device in 
accordance with the present invention. 

Description of the Invention 

Figure 1a illustrates a low temperature superconducting gradiometer 10 
in accordance with the present invention. Gradiometer 10 comprises a first 
field vector sensor in the form of superconducting wire loop 11. Gradiometer 10 
further comprises a second field vector sensor in the form of superconducting 
wire loop 12. Sensors 1 1 and 12 are positioned substantially coaxially so as to 
sense field vectors along axis 13. . Further, sensors 11 and 12 are connected 
in a differencing arrangement, whereby a signal produced at terminals 14 
stems from the difference between the magnetic field vectors sensed by 
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sensors 11 and 12. In accordance with the present invention, gradiometer 10 
is configured to rotate about axis 15, which is at a normal to the page as 
viewed in Figure 1a. Rotation of the gradiometer is indicated by direction 
arrows 16. 

5 Figure 1b illustrates a high temperature superconducting gradiometer 20 

in accordance with the present invention. Gradiometer 20 comprises a first 
field vector sensor in the form of a SQUID 21. Gradiometer 10 further 
comprises a second field vector sensor in the form of a SQUID 22. SQUIDs 21 
and 22 are mounted on a suitable magnetically transparent mount. Sensors 21 

io and 22 are positioned substantially coaxially so as to sense field vectors along 
axis 23. Further, sensors 21 and 22 are connected in a differencing 
arrangement (not shown), whereby a signal thus produced stems from the 
difference between the magnetic field vectors sensed by sensors 21 and 22. In 
accordance with the present invention, gradiometer 20 is configured to rotate 

15 about axis 25. 

Figure 2 illustrates a three-gradiometer magnetic field detection device 
30 in accordance with the present invention. The gradiometers are 
represented by cylinders 31, 32, 33. Cylinders 31, 32, 33 are respectively 
positioned on axes 41, 42, 43, and each of axes 41, 42 and 43 is not parallel 

20 with either of the other two axes. Each of gradiometers 31 , 32, 33 is rotatable 
about unique non-parallel axes in accordance with the present invention. 
Operation of such a device is set out in more detail in the following. 



25 



30 
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.This calculation models the sensitivity of the gradiometer instrument to the magnetic field, and. 
gradient tensor in the presence of 6 pairs of misalignment angles introduced , by; tilting, the;-. , 
SQUIDs through arbitrary angles. The effects of precession are not considered e^qplicitly but can 
probably be modelled to sufficient accuracy within the theory. It is demonstrated that the ■ 
misalignment errors lead to contamination of the gradient and B-field components by terms that . 
depend on the misalignments, but the effects of misalignment can be removed by the theory by 
introducing appropriate calibrations constants. A final solution is derived in terms of instrument 
coordinates and the transformations to geographic coordinates is also described. 

The resulting instrument provides a measurement of the true- value of the gradient , and field • 
components which is relatively robust with respect to changes in the calibration constants?. 
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I. OVERVIEW 

The .present embodiment comprises of a set of three rotating axial gradiometers or drum .grar\ 
diometers. In the 2-SQUID implementation the output voltages from a pair of SQUIDs operated'^ 
in flux-locked loops are differenced to provide information about the magnetic field density aid 
the first-order magnetic gradient tensor. If the SQUIDs are perfectly aligned with identical sen- 
sitivities then the difference voltage provides no information about the magnetic field density but 
this information is restored if the sensitivities of the SQUIDs are mismatched by a known amount. 
The difference voltage also provides a measurement of the first-order axial gradient provided that ^ 
the contributions from higher-order gradients are sufficiently small to be ignored. This is easily* 
seen by expanding the B-field as a Taylor series about the origin which is take to be centre point 
of the axial gradiometer. 



A. Taylor Series for functions of three variables. 

The Taylor series for a scalar function B{x,y,z) about the origin is 

In the case of a 3-dimensional vector field the previous equation applies to each of the vector 
components, so the a-component of B at the point r, B a (r) , and at the origin, B a (0) , are related- 
by 

3 a (r) = B a (0) + (x ± + y± + z§-^ B a + A (.£ + y£ 4- B„ 

3! \ dx y dy dz J 4\ \ dx dy dz J 

If the second- and higher-order gradients of B are.small then 2 may be approximated by 

-*■«>+ (*£ + »! + *JjK <3) 

provided that r is small. 
Let 

A 9 B 
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denote components of the the first-order gradients of B. Then eq. ?? becomes simply 

B a (r) == B a (0) + XQax + VQocy + ZQolz (4) 

or in coordinate independent form 

B(r) = B(0)+£r (5) 

where g is the gradient matrix and r is the displacement vector to the point at which the B field 
is measured. 

Using 5 if ri and r are two points the difference in the B-field at these points is 

B(n)-B(r )=B(0) + ^n-(B(0)+?r )=^(ri-r ; ) 

which is proportional to the product of the gradient matrix and the displacement between the 
points. 

II. INSTRUMENT AND DRUM COORDINATES. 

The drum gradiometers rotate about axes with orientations denoted by the vectors V "and" 
W in an umbrella configuration to be described in detail below. The axial gradiometers within the 
drums are located at the end points of U, V and W, and facilitate absolute value measurements 
of linear combinations of some of the components of the magnetic gradient tensor and magnetic 
field with respect to a drum coordinate system, C", which is fixed to the axis of rotation of- each* 
of the three drums. The systems of equations in drum coordinates are transformed to a common 
instrument coordinate system, C, and solved to yield absolute value measurements of the magnetic 
gradient tensor and magnetic field in instrument coordinates. 

For use in geophysical exploration the final measurements have to be expressed in terms of 
the geographic coordinate system, Co- This is easily done if the orientation of C with respect to 
Cq is known but in general this will require additional information from position and orientation, 
sensors. It is -possible!, . that, in some applications, the real- time measurements of the magnetic 
field provided by the ' instrument itself, might be used to provide information about the 

orientation of the instrument with respect to the magnetic field of the earth. Transformations of the 
final measurements between the instrument coordinates and geographic coordinates are described 
below in section VI. 
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A. Instrument coordinate system, C. 

The instrument coordinates system, C, is a right-handed Cartesian coordinate system defined 
as follows. The x-y plane is defined by the base of the instrument. The ^-direction is normal to 
the x-y plane and forms an axis of symmetry for the instrument in the sense that the three drum : 
gradiometers are distributed symmetrically about the z-axis. Unit vectors in the directions of the . 
x-, y- and 2-axes will be denoted i, j and k, respectively. 

As mentioned above the axial gradiometers are located at the end points of U, V and W and 
rotate about the axes defined by these vectors. The direction of U is the vector obtained by.' 
revolving a unit vector i with its base at the origin of C about the -y axis through, an angle A: , 
The directions of V and W are obtained by revolving XJ about the z-axis through angles of 120° 
and 240° respectively. U, V and W then form a type of umbrella configuration being disposed 
symmetrically about the *-axis and tilted above the x-y plane through an umbrella angle angle -Ai ' 
Hence 



U = (L u cos A, 0, L u sin A) (6) 
V = cos A, ^y~L v cos A, L v sin A^ (7) 

W= | — -jL^cos A, — —Lmcos X.Ln sin A I (8) 



(1 y/Z \ 

in which L u> L v and L w denote the distances of the axial gradiometers from the origin of C. 



1. Properties of U, V and W. 

• Note that these unit vectors form an orthogonal system only if 

U-V - V- W = W-U = 0 

i.e. if 

1 3 

7 - - cos 2A = 0 
4 4 



A = ^ arccos ~ = 0.615 48 rad = 35. 3° (9) 
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• In the instrument A = 45° so from equations 6, 7 and 8 

U=(coB^"lnj)=itt(M0M) ( 10 > 

V- (4 cos l^.cos Isinfj -^(-fc/S V*) ("> 

- (4 cos I>"f cos i> sin i) -^(-ivs V2) (12) 



w 



5. £frwt drum vectors 

It will prove useful to define the unit drum vectors, uo, vo and wo, that are normal to U, V 
and W respectively. Vector u 0 is the unit vector obtained by revolving a unit vector i with its 
base at the origin of C about the +y axis through an angle (90 — A) degrees. Vectors vo and wo 
are obtained by revolving u 0 about the 3-axis through angles of 120° and 240° respectively. In 
instrument coordinates 

uo = (sin A 0 -cos A ) (^ 3 ) 

(14)- 

: (15) 



wo 

N:B. 



vo==^-isinA -~~(sinA) — cosA^ 
= ^ — i sin A — ^-^(sinA) — cosA^ 



u 0 -U=(sinA 0 -cosA) • L u (cos A, 0, sin A) =0 

v 0 V= (-^sinA ^(sinA) - cos A ) ' L v cos A, ^ cos A, sin A^ =0 
wo • W= ^ S i n A (sin A) - cos A ) • L w cos A, cos A > sin A ) ===0 



as required. 



B. Drum coordinate systems C 



As mentioned earlier we introduce 3 drum coordinate systems, C", C" and C" associated withy 
the each of the three drums rotating about U, V and W respectively. Let x', y' and z' denote 
unit vectors in the directions of the x\y' and zf axes of the drum coordinate systems. C" is defined 
so that x' = uo, y' = j, z' = V/L u and the origins of C f v and C coincide. With these definitions 
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C f v is a right-handed coordinate system which would become equivalent to C if the umbrella angle 
was set to 90°. 

C" and C" are then defined by rotating C" about the z-axis through angles i20° and 240° so 
that in C", x" = v 0 and z" = V/L v while in C", x'" = w 0 and z"' = W/L w . 



1, Unit gradiometer vectors. 

The instantaneous direction of the rotating axial gradiometers will be defined, by unit gradiome- 
ter vectors u, v and w, which remain normal to U, V and W and undergo continuous rotation 
about the end points of U, V and W with angular frequencies u? u , u) v and oj w respectively. 6 Uy 
6 V and 0 W will denote respectively the angular displacements of the vectors u, v and w at time 
from the positions of u 0 , v 0 and w 0 defined above, so that 0 U = £J u t, 6 V = to v t and 0 W = w w t. In 
terms of the drum coordinates C' u% C f v and C f w the unit gradiometer vectors u', v' and w'rassume 
the same form 



u' = (cos 0 Ui sin 0 U) 0) 
v' = (cos 0 Vi sin 0 V , 0) 
w' = (cos 0 W) sin 0 Wi 0) 



(16) 
(18) 



C. Coordinate transformations between C and the C" systems. 

In general if P and P' denote the components of a vector in the frames G and C" respectively 
and if C is mapped to C by a positive (right-handed) rotation through an angle k about an axis 
k (k = x, y, z) then 

P' = ^(k)P ( 19 )- 

where R k (k) denotes the standard rotation matrices 

/l 0 

Rx (*0 = n cos k sin K I ( 20 ) 



0 cos k sin k 
y0 — sin/c cosk 

cosk 0 — sintt 



\ 



0 10 
\^ sin k 0 cos k J 



(2i) 
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R z (k) = 



V 



costt sin« 0 
— sinrc cos/c 0 
0 0 1 



\ 



(2?). 



Also if <7 and are the components of a Cartesian tensor in the frames C and C 9 respectively 
then 



r = Hfc(«)?i? r fc(K). 

For example a rotation about the -s-axis through angle k is 
1 9n 9i sfa\ 



9\ 9' 9'z 
\9l3 9' 3 933 J 
( 



cos k sin K 0 \ 
— sin k cosk 0 
0 0 1 J 



f 9X1 91 913^ 

91 9 9 3 
\913 9 3 933 J 



cos k — sin k 0 ^ 
sin k cos K 0 
0 0 1 / 



and hence 



9n = - (#11 + 9 ) + 1 (Pli —9 ) cos 2k + g\ sin 2k 
g' = - (#n + £ ) — - (gn - £ ) (cos 2k) - £1 sin 2k 
933 = 533 

£l = — - (gu — 0 ) (sin 2k) + pi cos 2k 
9* 3 ~ 0 3 cos k - P13 sin K 
#13 = 513COSK + P 3sin« 



(33)' 



(24), 



i. Thms/ormaiions C ^ C". 



When x' = uo,y' = j and z' = U/£ u , C is mapped onto C f by a rotation about y through an 
angle (n/2 - A) . Hence if B and B' denote the components of a vector in the frames C and C f 
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respectively then 



B' = 



/sin A 0 —cos A 

0 10 
y cos AO sin A 

( sin A 0 cos A 
0 10 



\ 



V 



— cos A 0 sin A 



\ 



B 



B 7 



(25) 



(26) : ; 



N.B. If A = tt/4 then 



f 1^/2 0 -V^ . .. ; ,. 

B'= 0 1 0 B (27)..; 

(Mo ^ j 

while if A = 7r/2 25 and 26 give that P' = P, i.e. the coordinate systems C and Cfa coincide, as 
required. Also, using equation 6, 

^ sin A 0 — cos A ^ 
U'= 0 1 .0 u 
^ cos AO sin A j 

f sin A 0 — cos A \ ( L u cos A ^ 



0 10 
cos A 0 sin A j 

0 

W 



^ L u sin A j 



as required. 

Using 26 the unit gradiometer vector 



u 



becomes 



u = 



sin A 0 cos A 



0 10 

cos A 0 sin A j 



( cos 0 U ^ 
sin 9 U 

V 0 J 

( q.osB u \ 
sin# n 

V 0 / 



( cos 9 U sin A 
sin 6 U 
— cos 0 U cos A 



(28) 
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in C* Observe 



cos 9 U sin A 



\ 



/ 



L u cos A 



\ 



u-TJ = S m0 u • 0 =0 

^ — cos 6 U cos X J \L U sin A j 

so u J_ TJ as required. 

a. Gradient Transformation. The transformations of the gradient tensor from C to C f follows* 
from the use of the transformation matrix in 25, with the general transformation equation, 23, i.e. 



^ 011 9i 013 ^ 

01 9* 9*3 
\9iz 9'z 033 / 



/ sin A 0 —cos A 
0 10 



\ 


^011 01 


013 ^ 


/ 


sin A 0 


cos A ^ 






01 0 


0 3 




0. 1 


0 : 


■ .'(29).; 


/ 


\ 013 0 3 


033 / 


\ 


— cos A 0 


sin A j 





011 = 1 (011 + 033 - 011 cos 2A -f #33 cos 2 A - 2gi3 sin 2A) 

0' *=0 

^ 033 = " (011 + 033 + 011 cos 2A - ^33 cos 2 A 4- 2^13 sin 2A) 
g[ = l V2{gx -93) 
013 = — - (2013 cos 2 A - gu sin 2 A + 033 sin 2 A) 
0 7 3 = 2 V2(0i +0 3) 

In the special case A = 7r/4 

011 = -011 - 013 + ^033 
9 1 =0 

033 = 2 0H + 013 + I 033 
01 =i V2(0l -0 3) 
013 = 2 0H - 2 033 
(7 ; 3 = ^V2(01 +0 3) 4 



(30). 



(31): 
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jg. Transformations C C n \ 



When x" = vo and z" = V/L v > O is mapped onto C" by a rotation about- z- through an angle 
27r/3 followed by a rotation about y through an angle (tt/2 — A) . 

^— -sinA -\/3(sinA) — cosA^ 



B" = 



B = 



-iy/3 -1 0 

^— -cos A --^(cosA) sin A y 
/ -A sin A -AV3 -±cosA \ 



V 



iv^CsinA) -i iV3(cosA) 
— cos A 0 sin A 



B" 



If A = tt/4 



B" = 



Using 33 the unit gradiometer vector 



B. 



^ cos 0 V \ 
sinO v 
0 



becomes 



( -^sinA -±V3 -^cosA \ 

^V3(sinA) -± IV^(cosA) 
^ —cos A 0 sin A y 



/ cos 9 V \ 
sin0 v 

v 0 



( — 1 sin A cos (sin0 v ) \ 

sin 0 V + ^ >/3 (sin A cos -ft,) 
^ — cos A cos 0 V 



in C. 

Observe that 



(32).; 



(33) 



(34). 



v V=0 
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confirming that vlVas required. Also using 6 

^-^sinA ^v^CsinA) -cosA^ 
V"^ -±v/3 - 1 0 V 

^ — - cos A -a/3 (cos A) sin A j 
''-isinA i%/3(sinA) -cosA^ 

^ — -cos A -V3(cosA) sin A J 



^ — -L v cos A ; 

^L v cos A 
^ L v sin A J 



= L V 



0 

v 1 ; 



as required. 

a. Gradient Transformation. The transformations of the gradient tensor from C to'C follows, 
from the use of the transformation matrix in 32, with the general transformation equation, 23, i.e. 



1 9ll 9\ 9l3 



\9l3 9 f 3 9h ) 



„lt Jfl „U 

9i 9 9 3 



V 



— icosA ^(cosA) sin A 



\ 


^ 511 91 513 ^ 


f -isinA 


_^ 




91 9 9 3 


^ (sin A) 


_1 


/ 


\ 913 9 3 933 J 


^ — cos A 


0 



sin A . 
(350 
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g'^ = l gn + | 5 + A333 - fan cos 2A - §3 cos 2A + ±333 cos 2A 
+i 5 i 3 sin2A - £31 \/3 + iffl \/3cos2A - ±g 3 \/3sin2A 

3" - iff + iffl V3 + V3 (iffl + iffllV3) 

333 = |ffu + |ff + 1 P33+iffiicos2A+|3 cos2A-i 53 3Cos2A 
-A313 sin 2A - |si V3 - |pi a/3 cos 2A + ±9 3^3 sin 2A 

"* 3i = 1 (cos A) 3 3 "(sin A) 01 + - (cos A) 3i3\/3 
+i(sinA)su\/3-i(sinA)3 \/3 

S J 3 = I ffl3 cos 2 A + ^311 sin 2A + §3 sin 2A - ^ sin 2A 
--9 3 \/3cos2A- 331 \/3sin2A 

^'3 = -i (cos A) si - 1 (sin A) 3 3 + \ (cos A) sil 
(cos A) 3 V3 - 1 (sin A) 3i3\/3 

In the special case A = Tr/4 

011 = iff" + i ffl3 + §3 + ^33 - i \/3~3l - 1 S 3V3 
3" = 3311 + \g + 1 3l V3 

333 = iffll - i Sl3 + Iff + 1 ff33~iV5ffl + V5ffS 
fl » =-|V23i +i\^S 3 + |\/6sil + iV63l3-iV63 

ffl3 = iffH + fff - A ff33-iV331 
g% = -\V2gi -iv^ps + iV^ffll-iV^fflS-iv^S 



(36) 



(37) 
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3. Transformations C C" 



If A = tt/4 



»//' . 



When x'" = w 0 and z'" = W/I- W) C is mapped onto C" by a rotation about * through an 
angle 47r/3 followed by a rotation about y through an angle (n/2 - A) . Hence 

^ -A sin A -A-s/3(sinA) -cosA^ 
1^/3 •-• 1 0 E 

— A cos A — iV3(cosA) sin A y 

-isinA A%/3 -icosA N 
-1^3 (sin A) - 1 -iv^(cosA) 
— cos A 0 sin A 



V 
/ 



B = 



V 



J 



B'" = 



Using 39 the unit gradiometer vector 



1^3 -1 Q 



cos N 
w"'= I sin^ 



(38); 



(39.). 



(40) 



becomes 



w = - 



-^sinA ±y/l -^cosA > \ 
^V3(sinA) - 1 -±y/Z(cos A) 
— cos A 0 sin A 



in C. 

Observe that 




( -± sin A cos 0 V + 1 V? (sin ^) ^ 
— i sin 0 V ~ - V3 (sin A cqs 0 v ) ' 
— cos A cos 0 v y. 



v 



(41) 



w W =0 
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confirming that wlWas required. Also using equation 8 





f- 


-sin A 


-iv^CsinA) 


— cos A ^ 








w"' = 








_1 


0 


W 








\- 


- cos A 


-i\/3(cosA) 


sin A j 












- sin A 


-±-\/3(sinA) 


— cos A ' 




/ 


—~L W cos A 


\ 








_1 


0 






-^L w cos A 






\- 


- cos A 


-^V^(cosA) 


sin A j 




V 


L w sin A 


/ 






(°) 
















0 


















K 1 ) 















as required. 

a. Gradient. Transformation. The transformations of the gradient tensor ftomGio-C^fo^ . 
lows from the use of the transformation matrix in 38, with the general transformation equation;, ; 
23, i.e. 



/ *jn jm ~w \ 

1 9\\ 9i 9i3 1 

9T 9 m 9% 

V«fi 9% S&J 



-isinA -^(sinA) -cos A 
& -I 0 

^ — 1 cos A — ^ (cos A) sin A j 



\ / 011 01 913^ 



91 9 9 3 
\913 9 3 933 } 



-^sinA ±y/3 -4cqsA 

_^( s inA) -1 ^4 (cos A 

— cos A 0 sin A 

(42) 
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fill - £su + h + ^33 - £sil cos 2A - |ff cos 2A + i<?33 cos 2A 
+i 513 sin2A+ \g\ y/3-\gi V3cos2A+±s 3 \/3sin2A 

533 = g5H + h + i 333 + |sil cos 2A + %g cos 2A - i<?33 cos 2A 
-i 5 l 3 sin2A+^5l V3 + |si V3cos2A - A s 3\/3sin2A 

** g>{' = i (cos A) g 3 - 1 (sin A) 51 - 1 (cos A) 513^ 
-^(sinA) ffU V3 + |(sinA)5 V5 

s ^ = 1 5l3 cos 2 A + Igu sin 2A + §s sin 2A 

-i$33sin2A + 1 g 3<y/3cos2A + 3SI v^su^A 

g% = -A (cos A) si - 1 (sin A) g 3 - \ (cos A) sil V3 
+^ (cos A) 0 -s/3 + A (sin A) 513^3 

In the special case A = ff/4 

511 = 1511 + 1 9U + h + A 533 + i\/3<71 + 1 V3<7 3 
ff"' = |511+i5 - 1 V351 

533 = &m-' i 513 + fs + 1 S33 + 3V^51 "^3 

g '» = -^V2si + J V2s 3 - iVesn - JVSsis + iVes 

S13=S511+|5 - 1 fl33+i>/Slw 

; - -iV2si - 3V^S 3 - I V69U + I n/6513 + iV6s 



9 3 



(43) 



(44) 



D. SQUID Configuration 



The SQUIDs composing the axial gradiometer rotating about U are located at the endpoints. 
of the vectors r = U ± pu, where p = d/2 and d is the total distance between the SQUIDs. In 
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drum coordinates 



/ o \ 



r' + = U' + pu' = 



r'_ = U' - pu' = 



p cos 6 U N 
/o sin 0 U 
0 

\ 



0 

\Lu) \ 
0 

\L U J \ 

Using eqs. 45 and 46 in equation 5 the magnetic field at the SQUIDs may be written 

( BL \ f g'xx 



p cos 6 U 
psin 6 U 
0 



p cos 0 U ^ 
p sin 8 U 

( —p cos 8 U y 
— psin 0 U 
L u 



B' + 04) = 



+ 



\ ( p cos 0 U ^ 



p sin # u 



<xy 9xz 
9yx a'yy 9y* 
\ 9zx 9yz 9zz J \ 
( B' x + Ly^ xz + p (sin d u ) g^y + P (c° s e u) 9x* ^ 

B' y + L u g' yz + p (cos 8 U ) g'^ + p (sin (?„) 
^ .B* + P (cos 0 U ) ^ + p (sin e u ) g' yz + Lug' zz J 



B'_ (r'_) - 



^ Pix 9'xy g'xz^ 
g'yx 9yy 9yz 



(b>\ 

B'y + 

\ B' z j \g'zx ff'zy g'zz ) 
( B' x + L u g' xz - p (sin 6 U ) g'xy-P (cos d u ) <4c ' 
B'y + L u g' yz - p (cos e u ) g^ - p (sin 9 U ) g'yy 



( — pcos# u 
— psin d u 



\ 



J 



V 



B' z - p (cos 6 U ) ^ xz - p (sin e u ) g' yz + L u o[ zz j 



in which 



(b>\ 

\B' Z J 



is the magnetic field at the origin of C. 



(45) 



(46);. 



(47)' 



(48). 



(49) 



E. SQUID outputs. 



In the instrument the gradiometer signal is produced by forming the sum 

U = U+ + U- 



(50) 
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where U+ and J7_ are the responses of the SQUIDs when operated'in a flux- locked loop. In response- 
to the magnetic fields in equations 47 and 48 U+ and U- may be written 

y + =si.B' + (r;)+a+ 

UL = S' -BL(rL)+a- . (52) : ' 

where the Si and S' denote the SQUID sensitivity vectors. These are normal to the SQUIDs with : < 
lengths proportional to the sensitivity of the SQUID system (V/ T) and the a are arbitrary offset . 
voltages which can change value upon each reset of the SQUID system electronics. In a perfectly 
aligned system Sj = Siu' and S' = — S u' but in any real system the SQUID Vectors will be tilted; . 
with respect to u' due to alignment errors. To characterize the alignment we introduce 2 pairs of 
alignment error angles. Let a\ (cr ) denote the error angle between the projection of the (S' )v 
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vector onto the x'-y* plane and the unit vectors u' (— u') , and t\ (r ) is the angle of inclination of 
S' x (S' ) above the x'-y f plane. Then 



si- 



S\ COS Tl COS 



S\ cos ti sin $i 
Si sin n 



, 0 A(6 u + <t ) 



\ 



J 



( S COST COS 6 

— S cost sin0 
S sinr 

Ifcri = <r ~ ri — r ==0 then = Siu' and S' = —5 u' as required. 
• To simplify these introduce 

Mi == Si cos n, iVi = Si sin n 
M ^ S cos r , iV = i9 sin r 



so that 53 and 54 become 



si- 



s' = 



/ Mi cos 
Ml sin (9j 

— M cos 6 
— M sin0 
iV • 



\ 



/ 



Using 47 and 56 in equation 51 gives 



Mi cos 0i 
Mi sin 0i 
ATl 



\ 



Y 



/ B' x + iu5^ + /=» ( sin 9 u) 9xy + /> ( cos 0 u) 9xx 
B' y + L u g' yz + p (cos 0„) g' xy + p (sin 0 U ) g'yy 
B' z + p (cos 0 U ) + p (sin 0 U ) ^ z + L u g' zz ) 

= (Mi cos 0i) (B' x + L u g' xz + p (sin 6 U ) g'^ + p (cos 0 U ) 
+ (Mi sin 6i) (B' y + L u g' yz + p (cos 6 U ) gt^ + p (sin 0 U ) g'yy) 
+ (JVi) (B' z + p (cos 0 U ) g' xz + p (sin 0 U ) g' yz + L u g' zz ) + a+ 

= B' x Mi cos 0i + L u g' xz Mi cos 0i + g' xx pMi (cos 0„) (cos 0i) + g'xypMi (sin 0 U ) (cos 0i) 
4- B' y Mi sin 0i + g' yz L u Mi sin 0 X 4- 5^pMi (sin 0 U ) (sin 0 X ) + p^Mi (cos 0 U ) (sin 0i) . 
+ B' z Ni + g' zz N x L u + g r xz N x p (cos 6 U ) + 9 yz XiP (sin 6 U ) + a+ 



(53). 



(54) 



(55) 



(56) 



(5.7) 



(58). 

(59) 
(60) 
(61). 



(62); 
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where we have used that = gyx, g xz = g zx and g yz = g^ from our knowledge of the properties 
of the magnetic gradient tensor. 

Similarly using 48 and 57 in equation 52 gives 

( -M cos6 \ / B' x + L u g' xz - P (sin9 u ) g'^-p (cos 9 u )g'xx^ 

-M sine ■ B' y + L u g' yz - P (cos 8 u )g' xy - p (sin 6 u )g'yy 
^ N J \B' z -p (cos 9 U ) g' xz - p (sin 9 U ) g' yz + Lvd„ ) 

= (-M cos 9 ) (B' x + Lug'** - P (sin 8 U ) g' xy - p (cos 9 U ) g'^) ' 
+ (-M sin 9 ) (By + L u g' yz - p (cos 0 M ) -g'^ - p (sin 0„) j^,) 
+ (JV ) (B; - /? (cos 9 U ) g' xz - p (sin 9 U ) g' yz + L u g' zz ) + <*- 



!7_ = 



+-0- 



(63) 

(64), 
(65) 
(66). 



= -B' X M cos 9 - Lug' xz M cos 9 +M pg xx (cos 9 U ) (cos 9 ) + M pg'^ (sin 9 U ) (cos 6 ) 
-M B' y (sin 0 ) - g' yz M L u (sin 0 ) + /> (sin 9 U ) (sin 0 ) + g'^M p (cos 0 U ) (sin 9 ) 

+ B^JV + g' zz N L u - g' xz N p (cos 9 U ) - g' yz N p (sin 9 U ) + a_ 



(67) 



WO 2004/015435 




PCT/AU2003/000999 



Substituting 62 and 67 into 50 gives the gradiometer output 

U = B' X M\ cos 0i + L u g xz M\ cos 0\ + g' xx pM\ (cos 0 U ) (cos 0\) 4- g^pMx (sin 0 U ) (cos 0i) 

- B X M cos 5 - L u g xz M cos 9 4- M pp^ (sin 0 U ) (cos J ) + M pp^ (cos 0 U ) (cos 6 .) 
+ J5^M X sin Ox 4 fi^^uAf i sin 0i + pJ^pMi (sin 0 tt ) (sin 9i) + g' xy pM 1 (cos 0 U ) (sin 0i) 

— M JBy (sin 0 ) - p^M L w (sin 0 ) 4- p^M p (cos 0 tt ) (sin 5 ) + g'^M p (sin 0 U ) (sinfl ) 
+ B^ATi + g r zz NiL u 4- flkr^ IP (cos 0 tt ) + p^JVi/o (sin 0 tt ) + a+ 

4- B*iV -f g' zz N L u - g' xz N p (cos 0 U ) - g f yz N p (sin 0 U ) + a_ "\ , 

= B' x Mi cos (0 U 4 cri) 4- X^p^Mi cos (0 U + <ri) + P^P^l ( cos ^) (cos (0 U ■+ cri)) 

+ PiyPMi (sin 0 U ) (cos (0 U + cri)) - B£M cos (0 U 4 a ) - L u g xz M cos (^ + a ) 

+ M pg'xy (sin fl tt ) (cos (0 U 4- cr )) + M />p^ (cos 0„) (cos (0 U + a )) 

4- B^Mi sin (0 U + <ri) 4- p^L u Mi sin (0 U 4 ai) + g'yyPMx (sin 0 tt ) (sin (0„ 4- <ri)") 

+ 9'xyPMi (cos 0 U ) (sin (0„ + cri)) — M By (sin (0 tt + <r )) 

- g' yz M L u (sin + a )) + p^M- p (cos 9 U ) (sin (0 tt 4 cr )) 

+ 4M p (sin 0 U ) (sin (0 U + cr )) + B^i + g zz NiL u 4 g f xz NxP (cos 0 U ) 

+ 9yzNlP (sin 0 U ) + a+ .+ B^iV 4- p^AT L u 

~ PfcoS *u) - g'yz N P ( Sin *u) + a - 

= B^Mi (cos cri cos 0 U — sin cr\ sin 0 U ) 4- L u p^Mi (cos ax cos 0 U — sin cri sin 0 U ) 

4- gt^pMi (cos cri cos 0 U — sin cri cos 0 U sin 0 U ) 4- gi^pMi (cos cri cos 0 U sin 0 U — sin <ri sin 0 U ) 

— B' X M (coscr cos 0 U — sin cr sin 9 U ) - L u g f xz M (cos a cos 0 U — sin cr sin0 u ) 

4- M pg xy (cos cr cos 0 U sin 0 U — sin cr sin 0 U ) 4 M pp^ (cos a cos 0 U — sin cr cos 0 U sin 6 U ) 

4- B^Mi (cos cri sin 9 U -4 sin cri cos 0 U ) 4- Py Z L u Mi (cos cri sin 9 U 4- sin <t\ cos 0 U ) 

4- g'yypMx (sin cri cos 0 U sin 0 U + cos cri sin 0 U ) + (cos cri cos 0 U sin 0 U + sin cri cos 0^) . 

— M (cos cr sin 0 U 4- sin cr cos 0 U ) - p^M L u (cos a sin 0 U 4- sin cr cos 0 U ) 

4- g' xy M p (cos cr cos 0 U sin 0 U + sin cr cos 0 U ) 4 p (sin cr cos 0 U sin 0 U 4* cos cr sin 0 U ) 

4- B^iNTi + g zz NiL u 4- ^ViP (cos 0 U ) 4- Py^iP (sin (9 U ) + a+ 
+ B^AT 4- p^iNT L u - g' xz N p (cos 0 U ) - g' yz N p (sin 0 U ) + 
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= B'xMi cos ax cos 0 U ~ B f x Mi sin <x\ sin 0 U 4 L u g xz M\ cos cri cos 9 U 

" &u9xz M i sin <*I sin 0 U 4- ^pMi cos a\ cos 0 U - g^pMi sin cri cos 0 W sin 0 U 

+ g'xypMi cos ai cos 9 U sin 0 U - ^pM x sin cri sin 0 W - jB^M coscr cos0 u 

4#iM sin cr sin 9 U — L u g' xz M cos a cos 9 U + L u g' xz M sin a sin0 u 

+ ^ Pd'xy cos ^ cos sin 9 U — M sin cr sin 0 U + M pg xx cos cr cos 0 W 

- M P9xx sin a c °s 9 U sin 0 U 4- B^Mi cos <x\ sin 0 U 4- B f y M\ sin cri cos 0 U 

+ g r yz L u M\ cos cri sin 6 U 4- p^Z^Mi sin cri cos 9 U 4 g^pMx sin cri cos 0<± sin 0^- 
+ 9vvP Ml cos a l sin &u + tfiy/'-Mi cos cri cos 0 U sin 0 U + g^pMx sin cri cos 0^ - 
-M By coscr sin 0 U - M sin cr cos# u -^ z Af L u coscr sin0 u 

— Pyjg-^ L u sin cr cos 0 U 4- g xy M p cos cr cos 9 U sin 0 U 4- #iy-^f P sin cr cos 0^ 
4p^ y Mpsincr cos 0 U sin 0^ 4 g'yyM«p cos cr sin 0 U 

4 B^iVi 4- g zz NiL u + c^Nip (cos 0 tt ) + g f yz N iP (sin 0 U ) 4- a+ 
4 B f z N -f L w - g f xz N p (cos 0„) - g f yz N p (sin 0 U ) + a_ 
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= B' x Mi cos 0 U cos <j\ - B f x M cos 6 U cos cr 
4- #£M sin a sin 6 U — -B^ Af i sin <ri sin 9 U 
4- ByM\ cos a\ sin 0 U — B^M cos a sin 0 U 
4- ByMi sin cri cos 0 U — B' y M sin cr cos 6 U 
4- I/tip^Mi cos cri cos 0 U — L u p^M cos cr cos 0 U 

— Lttp^Mi sin cri sin 0 U 4- Lug' xz M sin cr sin0 u 
4- c/y Z £ u Mi cos <ji sin 0 U — g' yz M L u cos cr sin 0 U 
4- ^i^uMi sin cri cos 6 U — g f yz M L u sin cr cos # u 

4- B' Z N! 4- B^TV 4- g'zzNxLu 4- tf^JV L u 4 a+ 4- a_ 
4- <^ x/ oMi cos a\ cos 0 U 4- M cos a cos 0 U 
4- g'yypMi cos cri sin 0 W 4 g'yyM p cos cr sin 0 U 
+ gvypM\ sin cri cos 0 U 4 sin cr cos # u 

— g'xypMi sin cri sin 6 U — M pg^ sin cr sin 0 U 

— g' xx pMi sin cri cos 0 W sin 0 U — M 1 pg^ sin cr cos 0 U sin 9 U 
4- g f yypM\ sin cri cos 0 U sin 0 K 4 pJ^M p sin cr cos 0 U sin 0 U 

+ Zg'xyPMi cos cri cos 0 U sin 9 U 4 2M ppiy cos cr cos 0 U sin 0 U 
+ 9xzN\P (cos 0 U ) - p^iV /o (cos 0 U ) 
+ g'yzNxp (sin 0 U ) - g f yz N p (sin 0 U ) 
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= (Mi cos a\ - M cos a ) B' x cos 0 U 
+ (M sin <r - Mi sin ai) sin 0 U 
4- (Mi cos &i — M cos cr ) sin 0 U 
4- (Mi sin — M sin cr ) B y cos 0 U 
4- X u (Mi cos a\ — M cos cr ) cos 0 U 

— L u (Mi sinai - M sincr ) #^sin0 u 
+ L u (Mi coscri — M coscr ) c/^sin0 u 

' 4- L u (Mi sin (Ti - M sin cr ) g' yz cos 0 U 
+ li»(JVi + ^)^ 
+ Si(JVi+JV) 
+ a 

4- /? (Mi cos 01 -f M cos cr ) ^ (cos 0 U ) 
4- /? (Mi cos cri + M cos cr ) g y (sin 0 U ) 
4- /o (Mi sin ai 4- M sin cr ) g' xy (cos 6 U ) 

— p (Mi sin ai 4~ M sin cr ) g'^ (sin 0 U ) 

— p (Mi sin c\ + M sine ) ^ (sin #u) (cos 0 U ) 
4- p (Mi sin <7i -h M sin cr ) g y (sin 6 U ) (cos 0 U ) 
4- 2/? (Mi cos cri -f M cos cr ) gxy (sin 0 U ) (cos 0 U ) 
+ «£*(cos0„) (N X ~N) 

+ P9'y Z (sin0 u )(N 1 -N) 

where a = a+ 4- a_. Expanding the trigonometric functions within each term and collecting 
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harmonics gives 

U= {L u g zz + B' z ) (N x + N) + a 
+ |p(Micosc7i + M coscr )(g' xx + g' yy ) 

+ (Mi cos cr\ - M cos a ) (B x + L u g' xz ) cos 0 U 
+ (Mi sin o*i — M sin cr ) + L u g yz ) cos 0 U 
+ p(JVi-iV)^cos^ u 

- (Mi sin a x - M sin a ) (J5^ 4- L u ^) sin 0 U 
4- (Mi cos g\ - M cos a) (B y + L u g' yz ) sin 0 U 
+ p(N 1 -N)g' yz sme u 

+ (^l cos at + M cos a- ) - p^) cos 20 u 
-f p (Mi sin cri + M sin a ) g xy cos 20 u 



+ 2^ ( Ml sinai + M sin <r ) (p^ - g xx ) sin20 u 
+ /? (Mi cos cr\ + M cos cr ) g xy sin 20 u , 



(68) 
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Hence 

U = \i u (L u g' zz + B' z )+a 

+ f « (g'xx + 9'yy) 

4- <5u (B x + L u g f xz ) cos 0 U 
+ Vu (By + L u g f yz ) cos 0 U 
+ v«Pas* cos 0 U 

- 4- L u g xz ) sin 0 U 

+ (By + I^u5^) sin 6 U 
+ ^^sin^ 

+ fu(s£xc-0^) cos 20 u 
+ 2 (Xu)9'xy cos 2G U 



+ Xu(g' yy -g'xx) sin20u 
+ 2 sin 20„. 



(69) 



where 



£ u = -ip (Mi cos a\ + M cos cr ) 
Xu — ~P (Mi sin 0*1 4- M sin cr ) 

^ (Mi cos a\—M cos cr ) 
Tfa = (Misincri — M sincr ) 
Mu=(JVi + i\0 

N.B. In the absence of alignment errors, i.e. a\ — a = 0, r\ = r =0 from equation 55 

Mi = 5i, iVi = 0 

M = S > N =0 
and from equations 70 Xu, *?u, and ^vanish, and 



(70). 



= ^(51 + 5) 1 
= (5l-5). J 



(71). 
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Finally if Si = S == S say, then <; u also vanishes and the only non-zero constant is . 

h-pS. (72). 

III. DATA EXTRACTION 

One method to extract the data is to take a series of measurements over a number of angular 
increments of the pick-up loop structure and solve the resulting system of equations algebraically. 
A series of voltage measurements is now taken at angles 0 U = — 7r/4, 0, 7r/4, ?r/2, 37r/4, tt, 57t/4,. 
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g' Xu + -PtiSaaMu. 

U ±U(G = tt/4) = a + B'zixu + giiCu - 9nXu + 2$i £« + ff' + 5' + iuff33Mu 
-iBfruVz + -B' Vu V2 + ^B'tfuVl + ±B' c«\/2 + VisWS + V 3 V «V^ 

U 4 Au(e = 3tt/4) = a + B&/x u + s' n ^ + ff'uXtt - 2pi + 9' U~ 9' Xu + L u g' 32 fJ. u 
"BirjuV^ - 1 B' T} U V2 ~ iB'^uVl + ±-B' c uV ^ - Vk^uV^ + V 3 *W2 
- 1 L u g' 13 r] u V2- i L u g , 3 r]uV2-^L u g , 13 <;uy/2+ 1 L u g' 3 <; u y/2 

Ue 'Au(0 = 5tt/4) = a 4- £ W + 0u?t* - 0uXi* + 2 0l f u + S>' £u + Xu + ^u^33^ . ' 

4- 1 Wi 3 t?u\/2 - ±JW 3 7? U V2 - li^0fc(uV5 - 1 L v ^ s ^iy/2 

Ux Au(6 = 0) = a + + B^u 4- + 2^ n ^ tt 4- si 3 v„ 4- 2g[ X u 

+L u g' 3 7] u + Lugitfu 4- L u g r 3S fj, u 

U z ±U(0 = tt/2) = a- Bi?7 u + £' c u 4- ~ Xu + s's^u + V £« 

— £ti<7i 3 77it + Lu9* zSu + ^u533Mu 

Us — U {0 = tt) = a — B l T] u — B[<; u 4- + 2^ tt - 013^ + 2 *i Xu 

Ur Au(6= 3tt/2) =a4 B^u - B'^ 4- B&u ~ 2g{ Xu - 9* *v u + 2^' 
4-£u0i 3 77u — ^lip 7 3^ + -ku033Mu 
Note that even with all of these misalignment errors considered the following L.C. vanishes if 
the SQUID stays in lock throughout the measurement series: 

Uq-Ui + U -l7 3 + £74-t/5 + %-^7=0 • / (74);. 

So this is still a good check to test the SQUID lock and the integrity of the data with respect- to/. 



3tt/2 to yield the following set of equations: 

Uo £U(e= -tt/4) = a + B&u + 9ii€u + 9iiXu-29i $ u + g* £ u - 
4- 1 J5 , 1 ?7u\/ 2 4- l B f <q u ^ 4- 1 B{<; u </2 - 1 B' <Tu\/2 + Vi3^u\/2 - 
4- i lWi3'7uV / 2 4- 1 Lu^ / 3 7 7uV2+ iiu^uV 2 - ^L^g' &*y/2 
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random pointing errors which might occur during the measurement sequence. 
A, Solving for gradient components. 

Assuming that Xu t Cu, \*>u are all known by calibration, the 8 equations in 73 contain 

10 unknown independent variables, namely B f x1 B f yi B* z , 

9xx> 9'yy, 9zz> 9xy> 9y Z > 9zx> °>\ sihoe voltages 
Uo„Uj are determined by measurement. The set of equations in 73 is clearly- uriderdetermined 
and additional data from the other two drums is necessary to determine a complete ; solution.' 
Nevertheless it is possible to solve for the off-diagonal gradient component, g f xyi and the difference , 
between the diagonal components, g 1 ^ — g'yy, independently of the measurements 'from the other 
two drums. This is done by forming the following linear combinations of eqs. 73 to efiminate the 
dependence on the SQUID offset voltage a and on the field components, 



Ui - U 3 + U 5 - U 7 « 4| u A' -f 8 X u9xy, 



(75) 
(76) 



where 



Then 



( 9 1 ^ 



A i A „' J 
— 9xx — 9yy- 



( XujUi -U3 + U5- U7) - U(Uo - U + Ua - Ue) \ 
8(£„ + xJ 

Xu(Uo -U +U4-U 6 ) + £ U (Ui -U3 + U5- U7) 
V " " 4(^ u + xJ J 



■ (77) 



provided 



(6u+Xu) ^0 



(78) 



which is always satisfied because both £ u and Xu are real and £ u > 0- 

The results from these additional gradiometers give solutions with exactly the same form in 
their respective coordinate systems, i.e. 

/ Xv(yi-Vz + V5-V7)-Zv(Vo-V + V4-V6) \ 
8(£*+Xv) 



V A"y 



XviYo-V +V4-V6)+$v(Vl-V3 + Vs-V7) 
4(^ + X«) 



(79) 
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and 



/ X W (^i — W3 4- W$ — W 7 ) -Z w (Wo-W +Wa- We) A 

8(^ + x w ) 



(80) 



X^O^O-T^ +^4-W 6 )+g w (Wi- W3 + W5-WV) 
V 4(*w + X«) 

Equations 77, 79 and 80 are now transformed to the instrument coordinate system and the. 

resulting system of equations may be solved to yield the gradient measurements. 



1. Transformations for the "proof of concept" instrument. 



In the case of the present embodiment A = ir/4, and the transformations for the xy\.xx and 
yy gradient components of the C", C" and C'" coordinate systems are those given by the; relevant' 
eauations in 36, 37 and 44, which are brought together here for convenience: 

g' x = Vafol -g 3) 

-\V2gi + lV2g 3 + lV6gn + \V6gi3 - §\/6ff 
-\V2gi + \-j2g 3 - gVosn - \V6gi3 + lV6g 



9'i 



9'{' 



9'n 



-011 — 913 + ~933 



9' =9 



g" = f 011 + \g + 1 gi V3 



9n 



(81).. 



i<m + hiz + Is + A S33 + iVsgi + 1 V3g 3 

s"' = |<7ll+30 -^Sl 
Using the last six of these equations we can derive the transformation equations for the diagonal' 
component differences A', A" and A"', i.e. 

&' = 9'u-9' 



= 2^11 -913-9 + -533 

1 A 1. 



(82) 



Aft ^// II 



5 1 1 3 /r 1 1 

= -g^ll + + -533 - ^01 v3 - -g 3 V3 + -5 13 

1a 1a 3 r- 1 /yl 
= --Ai + - A3 - -91 V3 - -5 3V3 + -513 



(83). 
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Ant _ jii mi 
= 9ii~ 9 

5 1 .1 3 y- 1 n 1 
= -gPil + + 2^33 + j9l v3 + -g 3V3 + -$13 

113 11 
= ~g Ai + -A 3 + -pi V3 + -<? 3 \/3 + - ffl3 (84) 

where 

A ±g - 533 > (85) 
A3 — 033 - 011 J 

Combining these with the first three equation of 81 yields the following set . of equations 

gi =iy^ 5l -iV2g 3 

9" =-|\/251 +\V2g 3+lV6gi3 + lV6Ai 

9i = -\V2gi + \y/2g 3 - \^giz - $>/6Ai I '. . . 

a' = -Ai — -a - 913 I , y ;; .;/v=-y%:. 

A" = -|Ai + iA 3 - fv/Isi - A V5<7 3 + l gi3 
A'" = -gAi + *A 3 + IVBgi + ^VSg 3 + 1 gi 3 4 

Although the corresponding matrix of coefficients, with respect to the variables 
gi .0 3. 013. Ai, A , A3, 







A V2 


0 


0 0 


0 \ 






f-0i \ 








-|V6 


0 


0 


5 3 




-*i 




IV2 


-\V2 




|\/6 0 


0 


513 




-ffi 




0 


0 


1 


_1 I 


0 


Ai 




"Ai 




|V3 


iVa 


_1 


i 0 


_ 1 


A 




-A' 


V 


-1V3 




_2 


ft 0 











is singular a solution can by obtained after observing that, by definition, 



Ai + A -f A3 = g n - g + g — 033 + £33 — 511 = 0 
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and use this equation to augment 86, i.e. 

Si = l V2gi - l <j2g 3 
g'{ = -\V2gi + \ V2g 3 + \V6gw + i\/6Ai 
g'{' = -\V5gi +l\/2ga- |V6si3 - |V6Ai 
A' = iA x -iA -513 

A" = -|A! + iA 3 - |\/3 5 i - i V3g 3 + i 5i3 
A" 7 = -|Ai + A A3 + |V3 ff i + ±V3g 3 + his 
0 = Ai + A +A 3 

The corresponding matrix form 



/-AV2 


AV2 


0 


.0 


0 


0 \ 


g 3 

913 

Ai 
A 




(~9{ \ 


IV2 








0 


0 






|V2 
0 


-\V2 
0 


1 


iV5 


0 
1 


0 
0 




~9\ 
-A' 




A^ 


_1 


1 

8 


0 


„ 1 




-A" 


-|V3 

V 0 


-A^ 
0 


_1 
0 


1 
8 

-1 


0 

-1 


-1/ 


\a 3 J 




-A'" 
I 0 / 



is transformed to a square matrix by multiplying on the left by the transpose of the coefficient 
matrix to give 



/ 33 
8 


3 


0 


0 


0 


0 ) 




(91 \ 


/ 


Vi - is? V2 - 1 A" V3 - Ipf v/2 + | A"'v^ > 


3 


9 


0 


0 


0 


0 




9 3 




-Vi V2+ Ja{ V2 - iA"V^+ i 5 f V2+;AA'"n/3 


0 


0 


9 
4 


1 

-4 


1 


1 




913 




-a' + i A" + * A'" + ^1 Ve - y{' V6 . 


0 


0 


1 


47 
3~ 


3 
I 


7 
8 




Ai 




1 a' - |a" - 1 A w + ve - ye V6 


0 


0 


1 


3 


5 
4 


1 




A 




-AA' 


V 0 


0 


1 


I 


1 








V 


A A" + -lA'" . / 



which is now non-singular. The inverse of the coefficient matrix is 





-1 


16 


0 


0 


0 


0 \ 




-« 


ft 


0 


0 


0 


0 




0 


0 


136 

5 


88 




6 

— 5 




0 


0 


88 


2Q4 




44 

-js 




0 


0 


64 
5 


11 
5 


436 


-% 


V 


0 


0 


-4 


_44 
75 


-ft 





WO 2004/015435 



41 



PCT/AU2003/000999 



and the solution of the system of equations, 86, is 

/ (A'" - A") + (2g[ - g'{ - g>{') \ 



(87) 



93 -^(A"'-A") + £ (9'{ -2ffi 

913 = JB (A" - 2A' + A'") + $ (rf - si") 

Ax & A' - £ A" - ^ A'" + £rf V6 - ^flT 

A -ii A' - * A* - * A'" - V5 + 

Note that at this point this solution has been obtained without using the tracelesshes's di *,the ■ 
gradient tensor and is written in terms of the differences between the. diagonal components. The. !; 
diagonal components, however, follow immediately from equations 85 and the traceless of tie 
gradient tensor, i.e. given that 

/ i _i n \ 



511 - 9 = Ai 
9 - 933 — A 
933 — 911 = A3 
911 + 9 + 533 = 0 



1-10 
0 1-1 
-1 0 1 
,111 



9 

\933 J 



A 
A3 

V 0 J 



(88) 



then 



( 9 ll\ 




^-A 3 + Ai^ 




9 


_ 1 
~~ 3 


-Ai+ A 


(89) 


\S33 J 




^ -A + A 3 J 





Substituting for the component differences from the last three equations in 87 gives 



(90) 



• ( -35 (-A' + 8A" + 8A'") + Wl - g'{') \ 
-3g (11 A + 2A" + 2A'") + ( P f - g'l ) 
§ (A' + A" + A"') . 

The final solution for the gradient components follow after substituting for the A's and.s's.: 
from equations 77, 79 and 80 into equation 90 and the first three equations of 87. The resulting' 



( 9 n\ 

9 

\933 J 



\ 



J 
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equations are brought together here for convenience, 

m = - u (-A' + 8A" + 8 A'") + ( 5 ? - *f ) 

9 =- M (HA' + 2A" + 2A'") + -£ (gl' - g'{ ) 

533 = i(A' + A" + A'") 

91 =^(A'"-A") + if (2si -<tf -Pi") 

P 3 = -# (A'" - A") 4- £ (p? - 2 5 i + pf ) 

513 = T5 (A" - 2A' + A'") + 4 (g'{ - g'i') 

in which 



, xuO/i -U3 + U5- Ur) -t u (Uo-U + U* - Ue) 
91 " 8(£ u + x,i) 



„ X„ (Vj - V 3 + V5 - W) - €v (Vb - y + V* - Ve) (91) 
51 " 8(^ + x«) 



P?' = 



V — 



A" = 



8 (t w + Xw) 

Xu (Up - U ± Ua - Ue) + tu (Ui - Ua + Us - U7) 
4(f„ + X u ) 

Xv(Vo-V +V4-Va) + ev(Vi-V3 + VB-V7) 
4«« + Xw) 



„, xw ( Wo - W + W 4 - W 6 ) + gw (Wi - W 3 + W 5 - Wg) 



a'" 

4(*v, + X«.) 



and 

^ u ^: ip (Mi cos <ri + M cos cr ) 
Xu — -P {Ml sin cr! + M sin cr ) 
Ml = 5i cos n, Ni = Si sin n 
M =S cost , TV =5 sinr 
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2. Special case: "Proof of concept instrument 11 with zero alignment errors. 
If all of the alignment errors vanish then 

Xu = 0 

and equation. 91 reduces to 

9X1 = -4s (-A' + 8A" + 8A'") + (g'{ - pf ) 

9 =- 3 g(llA' + 2A" + 2A'") + i^(fff "Pi) 
533 = | (A' + A" + A'") 

91 = (A'" - A") + (2 9 [ - g'{ - g>{>) 
9 3 = -# (A'" - A") + £ {g'{ - 2g> x + g>{>) 
913 = is (A" - 2 A' + A'") + 4 (g'{ - g'{> ) 

and (9?) 

, -(Uq-U +U4-U 6 ) • ... 

Pl " 8^ 

» -<Yo-V +V4- Ve) 

91 = 55 

w - (Wo - W +W A - W 6 ) 
91 = 8& 

V (Ul -U3 + U5- ui) 

4^ 
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IV. SOLVING FOR FIELD COMPONENTS, 

If the gradiometer is intentionally unbalanced, i.e. Si ^ S , then it is capable of returning 
information about the magnetic field in addition to the gradient. As for the gradient case considered 
above we form linear combinations of 73 to both eliminate the dependence on the SQUID offset 
voltage, a, and to reduce the number of gradient terms. To eliminate the SQUID offset we start 
with the differences between voltages corresponding to opposing positions, and then take the sums: 
and differences of these pairs, i.e. 

(Uq - 04) + (U - U 6 ) = 2rj u B' V2 + 2<; U B[ V2 + 2v u g f 13 V2 + 2L u r} u g f 3 V2 + 2L u <:u9xz^ 

(Uq - U4) - (U - Ue) = 2?7 u J5i\/2 - 2<; U B' V2 - 2^3^ + 2L u r)u9i3 ^ - 22^0' 3V2 . ; ■ * 

(U 1 -U 5 ) + {Ua-U 7 ) = 2B' Vu -2B' 1 r! u + 2B , l <; u + 2B'<; u 

+ 2g' 13 v u 4- 2p' 3 t; u — 2L u g' 13 r} u + 2L u g' 3 rj u 4- 2X^13^ 4- 2L u g' 3 q u 

(Ux-Us)-^- U 7 ) = 2^7^ + 2£' 77 U + 2£^ - 25' 

+ 2g' 13 v u — 2g f 3 v u + 2L lt ^ 3 77 u -f- 2L u g r 3 7] u + 2Ir u pi 3 Cu — 2L u g r 3 <; u 

and then simplify further by forming the following L.G. of these equations: 
Ui + Us-Us-Ur Uo + U - L/4 - i7 6 



2 2-n/2 
^1 - Us - Us + Cfr Uo-U -U 4 + U 6 



= CuB' - r) u Bi + vup'a - -Lu^ia + LuSug' z (93) 
= Tj u B / + <; U B[ + v u g f 13 + L u r}ug f 3 + Lrfu&a ' (94) 



2 2^2 
As for the gradient calculation previously the results from the other two drums have exactly the 
same form in their respective coordinate systems, i.e. ■ ' • 

Vi + V Z - V 5 - V 7 Vo + V - v 4 -v 6 



2 2y/2 

Vi - V z - V5 + V 7 V 0 - V - V4 + Ve 



= <; V B" - 7) V B'{ + v v g"z ~ L vVv g f { 3 + L v <^" 3 ' (95) 
= + <; V B'{ + ^'3 + W3 + (9 6 ) 



2 2V2 



2y/2 

+ ^iu0 7 3 LwVw9iz + L w <; w g m 3 (97) 
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W t -W 3 -W 5 + W 7 Wo-W - W4 + Wis _,,„ , D ,„ 
2 iVf = 7? ^ + ^ 1 

4- Vu,/ 3 + L w r) w g'" 3 + L w <; w g f { 3 

Equations 93 - ?? may be simplified by making the following substitutions: 

a £ ^-7= S 3 (^"Cw + v u ) + L u g' l3 r} u 

u 2 2J2 913 * " ^ ~ u9 3?7w 

2 2V2 L v9"zVv - tf zVv - L v g" 3 <; v 

K 2 ^-7= 9i3 v v - L v g" z r} v - Z-„p^c„ 

Tr t A W± + W * ~ ^ 5 ~ Wp + W - W 4 -W& 



-(98), 



2V2 



4- L w g 13 r) w — 3 3^10 — i>tu5 3^10 . 



(99) 



w 2 2vf S 3 U «> - L w g'" 3 T] w - L w g'{ 3 <; w 

Observe that all of the terms on the right hand sides of these equations are known i.e. the U, V 
and W voltages are returned by instrument, the gradient components are known in the instrument 
frame and can be easily found using coordinate transformations, and the remaining terms are 
known by calibration. With these substitutions 93 - ?? become 

tf + = «J3' - 77 u Si 

17* = t)uB' + CuBi 

Vt = s v B" - tj v B'{ 

V* = VvB" + s v B'{ 
W^ = s w B'" - T) W B'{' 
W* = r) w B'" + wBi" j 

Equations 100, corresponding to measurements in the C, C" and C" coordinate systems, may 
be solved. by transforming to the common instrument coordinate system using the coordinate-, 
transformations for field components derived earlier, but the resulting sets of equations rapidly 



(100) 
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become intractible. It is simpler to first solve these pairwise in their respective coordinate systems, 




(101) 



B'{ \ = 1 { V^ v - V* Vv 



. (102) : 



( 



(103) 



and then carry out the transformation to the instrument system. As derived above the relevant, 
transformation equations for the field components are 



B' = 



sin A 0 —cos A 



B'" = 



0 10 

cos AO sin A 

B'x = - (cos A) B 3 + (sin A) Bi 
B' — B 

B 3 = (cos A) Bi 4- (sin A) B 3 

''-A sin A lv/3(sinA) -cosA^ 

.1^3 0B 
^— -cos.A A\/3(cosA) sin A j 

' B'{ = - (cosA)B 3 - 1 (sinA) B x + * (sinA)B V3 
B" = -^B - ABin/3 

J3g = -1 (cos A) Bi 4- (sin A) B 3 + 1 (cos A) B V3 

/-A sin A —1^3 (sin A) -cosA^ 

1^3 -1 0B 

y^— -cos A — -%/3(cosA) sin A J 
' B'{' = - (cos A) B 3 - 1 (sin A) Bi - 1 (sin A)B s/% 
< B'" = — -B + iBiVS 

^ Bg' =-i (cos A) Bi + (sin A) B 3 - i (cos A)B \/3 
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and if A = 7r/4 



B'i= V2(£i-B 3 ) 
B' = B 

B'{ = \V2(-Bi+V3B -2B 3 ) 
B" = -i(BiV3 + S ) 
Bf = |v^(-Bi- V3B -2B 3 ) 
B"' = 1 (V^Bi - B ) 

and equations 101 - 103 become 



= ^(-Bi + \/3B -2B 3 ) 
= -i(BiV3 + B ) 



Vv + <*> 

y*»7 V + y*^ 



= |V5 (-Bx - V3B - 2S 3 ) 



(1.04) 



(105) 



Vw + ^ 

A least-squares best-fit solution may be found by first rewriting this over deter rained system of 
equations in matrix form in terms of the unknowns B\>B , £3 to give 



0 









( U^u-U^u \ 








Jlu + «u. 
U*T} U + U^u 


0 
0 


/b^ 

B 




Jlu + % 

vhv - v^vv 
vt Vv + vk v 




\Bz) 




Wk w - W^ w 


0 / 






.Vw ■ ">to , 

ty*??*, + 
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and then multiplying on the left by the transpose of the coefficient matrix to yield 

i 1 v k u ^ ~ u Vu - ~V2 - ha _ 1 * w * w ~ w ^ 
2 nu + Su 4 Vv + <^ 4 7} w + <^, 



1 ^ Whw-W^r} w _ lVhjv + V^ _ 1 W*r, w + W^ w 
4 T? U) + ? U , 2 7?„ + Sv 2 ^uj + Su, 
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We bring the final forms of the solutions here for convenience 

Uku ~ U*Vu V2 Vkv - V^ijy 2s/3 (V^v 4- VU V ) 
1 " 9 77 u -h C u 9 r7 v + C v 9 T7„ + S v 

_ V2 ^^V + 2>/3 (W^ffr, + g^tg) (1Q7) 

• iJw + i 9 Vw + S w . ' 

-D = — ; h 



9 ^4-^ 9 Vv + Sv 9 + 

9 + Cu, 9 Vw^^w 



where 



(108) 



^ A u k± u 1 -u 1 -u 1 _ Uo±u_^^ + Lug , 3Vu _j 3Vu _ Lugl ^ 

ut A u^-u^u, _ Uo-U-U 4 + Ue _ g , zvu _ LW3Vu _ WiaCu 

V T = — -7= h Lv9uVv ~ 9 3 v v — L v9 3«u 

w A V, - a - V 5 ± V 7 _ Vo - V-V< ± Ve _ ghvv _ Lvg „ 3Vv _ L ^ • 

^ ^ W 1 + W 3 -W 5 -W 7 _ Wo + W-W<-W 9 + _ ^ _ 

+ = 2^/2 3 "~ 3 iu, 5l3^tu 

. 513 = A 0U - 1 933, 9"z = ~3\/2pi - Jv^fl 3 + § V^il - {^913 ~ . 

g' 3 = ±V2(gx +fl 3 ), 013 = gSll + %9 - X 933 + \V^9l 

9i3 - Iffli + fs - i P33 - I V3ffl , 9" 3 = -\V*9l ~ \V2g 3 ~ 5 V6ffll + |\/6ffi3 +. \y/Bg . 
Formally the scaling coefficients, v u , ? u and r) u axe denned by equations 55 and .70 but would 
normally be determined by calibration. 
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EXAMPLES 



Precession about the 2-axis. 



Say a\ = a = 0, Si = S = S. If the gradiometer precesses about the the axis of rotation by- 
virtue of an error angle of cj> say, then the theory developed above suffices to model this by letting^ 
n = — r = <f>. 1 . Using equations ?? 



£ u — pQ£ cos 4> 
u u = 2pQS sin 



a u = (2Q S p cos (£) 1 

C/* = C/i — L/3 -4- C/s — t/7 
Mi^Scos<£, Nx^Ssuk/) 
M ± S cos cf>, N ^ -Ssm<j> 



(iiQ) 



VI. APPENDIX 1: GEOGRAPHIC COORDINATE SYSTEM, C G . 

The geographical coordinates system, C(?, is defined as follows: 
XG-axis: The a;c?-axis is along the traverse in the direction of increasing station number ( "grid north" ) . . 

2/Q-axis: The yo-axis is at 90° to the rcc?-axis, parallel to the ground and points right as you look" 
along the traverse ("grid east"). 

ZG-axis: The z<3-axis is vertically down (pointing at the centre of mass of the earth). 

This is a right-handed Cartesian coordinate system in terms of which a unit vector with azimuth. 
A (the angle measured clockwise from grid north) and inclination I (angle measured positive 
downward from the horizontal) is 



unit vector with azimuth A and inclination I = (cos A cos/, sin A cos I, sin I) . 



(in) 



1 If the gradiometer drum is tilted through and an angle </> in the x-z plane say, the pick-up loop vectors are also 
tilted through the same angle and the positions of the pick-up loops are altered. The change in position has a 
negligable effect on measurements of gradients produced by distance sources. Since by definition r is the angle' of 
the sensitivity vector above the x-y plane, a tilt though <t> implies that n = <f> and r = —<f>. 
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1. Transformations between C and Cg- 

The transformation between C and Cg assumes a particularly simple form if the instrument is 
oriented with its base horizontal and its x-axis aligned due north. In this case C is mapped to Cq 
by a positive rotation through an angle it radians about the x-axis and the transformations of the 
field and tensor and data become simply 



B G = RttB 
9G = R^gRn 1 



where 



fx o 

0 -1 
\0 0 



0 \ 
0 

- 1 / 



(112) 
(113) 



.(114) ' 



Hence 



and 



9G 1 9G 







[ Bl \ 








b g 




—B 






(115) 


\B G 3 J 




\- B *) 








9G13 ^ 




( 911 ' 


-91 


-913 ^ 




9G 3 




-9 1 


9 


9 3 


(116) 


9GZZ / 




^ -931 


93 


533 / 





More generally if the instrument has azimuth A and inclination I (in geographic coordinates)' 
then C is mapped to Cg by first rotating about y through an angle of —I then by a rotation about 
the resulting z through an angle A followed lastly by a rotation about the resulting x of 7r: Hence 



B G =R X (n)R z {A)Ry{-I)B 



(H7) 



f Bgi^ 
B G 

\BG3 J 



( 1 0 0 
0 cos 7r sin 7T 



0 — Sin7T C0S7T 



N ( cos A sin A 0^ 
— sin A cos A 0 
0 0 1 



cos A cos/ sin A cos A sin 



sin A cos J —cos A sin A sin J 
sin J 0 — cos / 



/ 

(b x \ 

B 

\B 3 J 



cos(-J) 0 -sin(-J) 



\ 



old 

^ sin (— J) 0 cos (— I) J. 



f Bl \ 

B 
\BiJ 



(118) 
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and 



/ 9G11 9GX 9G13 ^ 

9G1 9G 9G 3 
\ 9GX3 9G 3 9G33 / 

/ cos A cos I sin A cos -4 sin I ^ 

sin A cos J —cos A sin A sin/ 
^ sin/ 0 -cos I j 

In the special case A = J = 0 we get so that 

/ Bgi\ 

Bg 
\Bgz ) 



( 9X1 91 913^ 

91 9 9 3 
\913 9 3 #33 J 



-B 
\-BzJ 



( cos A cos/ sin A cos/ sin/ \ 

sin A — cos A . 0 
y cos A sin / sin A sin / — cos / J 



( 9GU 9G1 9G13 ^ 

9GX 9G 9G 3 
y 9GX3 9G 3 9G33 J 



-9X3 ^ 



9XX -91 
~9X 9 9 3 
\ —913 9 3 933 J 



as above. 
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The approximation used in obtaining equation (3) above, which is 
significantly simplified as compared to equation (2), involves truncating the 
Taylor series expansion of the magnetic field as a function of position to include 
the effects of the B-field and the first-order gradients only. This simplifying 

5 approximation is desirable because the Taylor series expansion has an infinite 
number of terms. The approximation is justified because the magnitudes of the 
higher-order terms decreases as the order of the terms increases, for three 
reasons. The higher-order terms in the Taylor series are scaled by prefactors 
which may be written as the product of two factors. Firstly, one of these factors 

to is a scaling term which decreases with increasing order. The second factor is 
proportional to 1/r", where r is half the distance of separation between the field 
detectors, and n is the order of the term. Since the separation between the 
field sensors is small, raising this to an increasing power rapidly decreases the 
magnitude of this factor. The third reason the approximation is justified is that 

15 higher-order gradients fall off faster with distance from the field source than 
lower-order gradients, so the signal at the gradiometer is dominated by the field 
and the first-order gradient components. 

The effects of misalignments of the field detectors on the output voltages 
from the instrument are also scaled by these higher-order terms, but these 

20 effects are small and similarly decrease with increasing order. 

Referring to equations (6)-(8), it is noted that the value of X may be any 
value in the range of 0° to 90°, inclusive. As noted in equation (9), when X is 
approximately 35.3°, an orthogonal system is provided. In the presently 
described embodiment, X is 45°. 

25 It will be appreciated by persons skilled in the art that numerous 

variations and/or modifications may be made to the invention as shown in the 
specific embodiments without departing from the spirit or scope of the invention 
as broadly described. In particular, the present invention may be applied in 
conjunction with the invention set out in the co-pending PCT application 

30 claiming priority from Australian Provisional Patent Application No. 2002950624 
and filed simultaneously with the present application, the contents of which are 
herein incorporated by reference. The present embodiments are, therefore, to 
be considered in all respects as illustrative and not restrictive. 
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CLAIMS 

1 . A method of field detection comprising: 

providing a gradiometer in a field, the gradiometer having at least first 
and second field vector sensors connected in a differencing arrangement; and 
5 controllably altering a position of the at least first and second field vector 

sensors relative to the field during operation of the gradiometer. 

2. The method of claim 1 wherein the at least first and second field vector 
sensors of the gradiometer are substantially axially aligned, such that the 
gradiometer is an axial gradiometer. 

io 3. The method of claim 1 or claim 2 wherein the step of controllably altering 
the position of the at least first and second field vector sensors is performed by 
rotating the at least first and second field vector sensors about an axis of 
rotation. 

4. The method of claim 3 wherein the step of controllably altering the 
15 position of the at least first and second field vector sensors is performed by 

rotating the at least first and second field vector sensors continuously during 
operation of the gradiometer. 

5. The method of claim 4, further comprising the step of providing a motor 
for driving the gradiometer at a substantially constant revolutional velocity 

20 about the axis of rotation. 

6. The method of claim 3 wherein the step of controllably altering the 
position of the at least first and second field vector sensors is performed by 
rotating the at least first and second field vectors piecewise about an axis of 
rotation. 

25 7. The method of claim 3 when dependent on claim 2 or any one of claims 
4 to 6 when dependent on claims 2 and 3, wherein the axis of rotation is 
positioned substantially perpendicular to the axial alignment of the first and 
second field vector sensors. 

8. The method of claim 7 wherein the axis of rotation is positioned 
30 substantially between the first and second field vector sensors. 

9. The method of claim 7 or claim 8, wherein the axis of rotation is 
positioned substantially equidistant from the first and second field vector 
sensors. 

10. The method of any one of the preceding claims wherein the gradiometer 
35 is an axial magnetic gradiometer. 
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11. The method of claim 10 wherein the field vector sensors of the axial 
magnetic gradiometer are SQUIDs. 

12. The method of claim 10 wherein the field vector sensors of the axial 
magnetic gradiometer are flux gates. 

5 13. The method of claim 10 wherein the field vector sensors of the axial 
magnetic gradiometer are superconducting pick-up loops. 

14. The method of any one of claims 10 to 13 wherein the sensitivity vectors 
of the field vector sensors lie substantially in a nominal x-y plane. 

15. The method of claim 14 wherein the axial magnetic gradiometer is 
10 rotated about a nominal z-axis perpendicular to the x-y plane. 

16. The method of any one of claims 10 to 15 further comprising the step of 
retrieving information relating to the magnetic field from the axial magnetic 
gradiometer as its position is controllably altered. 

17. The method of claim 16 wherein the information relating to the magnetic 
15 field comprises the gxy component and linear combinations of the vaiues of the 

gxx and gyy components of a field gradient tensor, in addition to information 
about the B x and B y field components. 

18. The method of claim 16 or claim 17 wherein the step of retrieving 
information relating to the magnetic field further comprises transforming a 

20 received signal from the field vector sensors into the Fourier domain in order to 
retrieve the information relating to the magnetic field. 

19. The method of any one of the preceding claims further comprising the 
steps of: 

providing second and third gradiometers, each of the gradiometers 
25 having at least first and second axially aligned field vector sensors connected in 
a differencing arrangement; 

positioning the gradiometers such that an axis of each gradiometer is not 
parallel to an axis of any other of the gradiometers; and 

controllably altering a position of the at least first and second field vector 
30 sensors of each gradiometer relative to the field during operation of that., 
gradiometer. 

20. A field detection device comprising: 

a gradiometer, the gradiometer having at least first and second fieid 
vector sensors connected in a differencing arrangement; and 
35 means for controllably altering the position of the at least first and second 

field vector sensors relative to a field during operation of the gradiometer. 
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21 . The device of claim 20, wherein the at least first and second field vector 
sensors of the gradiometer are substantially axially aligned, such that the 
gradiometer is an axial gradiometer 

22. The device of claim 20 or claim 21, wherein the means for controllably 
5 altering the position of the at least first and second field vector sensors 

comprises means for rotating the at least first and second field vector sensors 
about an axis of rotation. 

23. The device of claim 22 wherein the means for rotating the at least first 
and second field vector sensors about the axis of rotation is operable to rotate 

io the at least first and second field vector sensors continuously during operation 
of the gradiometer. 

24. The device of claim 23 wherein the means for rotating is a motor 
operable to drive the gradiometer at a substantially constant revolutional. 
velocity about the axis of rotation. 

is 25. The device of claim 22, wherein the means for rotating the at least first 
and second field vector sensors about the axis of rotation is operable to rotate 
the gradiometer piecewise about the axis of rotation. 

26. The device of claim 22 when dependent on claim 21 or any one of claims 
23 to 25 when dependent on claims 21 and 22, wherein the axis of rotation is 

20 positioned substantially perpendicular to the co-axial first and second field 
vectors. 

27. The device of claim 26 wherein the axis of rotation is positioned 
substantially between the first and second field vector sensors. 

28. The device of claim 26 or claim 27 wherein the axis of rotation is 
25 positioned substantially equidistant from the first and second field vector 

sensors. 

29. The device of any one of claims 20 to 18, wherein the gradiometer is an 
axial gradiometer. 

30. The device of claim 29, wherein the field vector sensors of the magnetic 
30 axial gradiometer are SQUIDs. 

31 . The device of claim 29, wherein the field vector sensors of the magnetic 
axial gradiometer are flux gates. 

32. The device of claim 29, wherein the field vector sensors of the magnetic 
axial gradiometer are superconducting pick-up loops. . . 

35 33. The device of any one of claims 29 to 32 wherein the sensitivity vectors 
of the field vector sensors lie substantially in a nominal x-y plane. 
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34. The device of claim 33 wherein the means for controllably altering the 
position of the at least first and second field vector sensors comprises means 
for rotating the at least first and second field vector sensors about a nominal z- 
axis perpendicular to the x-y plane. 
5 35. The device of any one of claims 29 to 34 further comprising means for 
retrieving information relating to the magnetic field from the axial magnetic, 
gradiometer as its position is controllably altered. 

36. The device of claim 35 wherein the information relating to the magnetic 
field comprises the gxy component and linear combinations of the values of the 

io gxx and gyy components of a field gradient tensor, in addition to information 
about the B x and B y field components. 

37. The device of claim 35 or claim 36 further comprising means for 
transforming a received signal from the field vector sensors into the Fourier 
domain in order to retrieve the information relating to the magnetic field. 

15 38. The device of any one of claims 20 to 37, further comprising: 

second and third gradiometers, each of the gradiometers having at least 

first and second field vector sensors connected in a differencing arrangement, 

the gradiometers being positioned such that an axis of each gradiometer is not 

parallel to an axis of any other of the gradiometers; and 
20 means for controllably altering a position of the at least first and second 

field vector sensors of each gradiometer relative to the field during operation of 

that gradiometer. 

39. A method of obtaining a complete magnetic gradient tensor of a 

magnetic field, the method comprising: 
25 providing at least three axial gradiometers in the magnetic field such that 

an axis of each axial gradiometer is not parallel to an axis of any other one of 

the at least three gradiometers; and 

controllably altering a position of each of the at least three gradiometers 

relative to the magnetic field during operation of the gradiometer. 
30 40. The method of claim 39 wherein the step of controllably altering 

comprises rotating each of the at least three axial gradiometers about a 

respective axis of rotation of that gradiometer. 

41 . The method of claim 40 wherein the axis of rotation of each gradiometer 
is not parallel to an axis of rotation of any other one of the at least three 
35 gradiometers. 
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42. The method of claim 41 or claim 42 wherein the at least three 
gradiometers are rotated at a substantially constant velocity about their 
respective axes of rotation. 

43. The method of claim 42 further comprising the step of transforming 
5 obtained data to the Fourier domain. 

44. The method of claim 43 further comprising the step of distinguishing field 
gradient information from field information in the Fourier domain. 

45. The method of claim 43 or claim 44 further comprising distinguishing 
information about the gxy component of the gradient tensor from information 

10 due to the diagonal gradient components, even at the same frequency. 

46. The method of any one of claims 43 to 45 wherein the Fourier transform 
is performed by a fast Fourier transform algorithm. 

47. The method of any one of claims 43 to 45 wherein the Fourier transform 
is performed by a simple discrete implementation via a coarse sampling of a 

15 discrete number of data points per cycle of the rotation. 

48. The method of claim 47 wherein eight data points are sampled. 

49. The method of claim 47 or claim 48 wherein the data points are 
positioned equally around the 360 degree rotation. 

50. The method of any one of claims 42 to 49 further comprising the step of 
20 applying a least-squares best-fit solution to obtained data. 

51. The method of any one of claims 42 to 50 wherein DC offsets are 
determined and monitored to provide information about the operating 
conditions of the gradiometers. 

52. The method of claim 51 wherein the DC offsets comprise one or both of: 
25 low frequency drift in at least one field vector sensor of the at least three 

gradiometers; and the fixed offset of at least one field vector sensor of the at 
least three gradiometers. 

53. The method of any one of claims 42 to 52, further comprising the step of 
isolating the at least three gradiometers from physical vibration occurring in the 

30 vicinity of their respective fundamental and first harmonic frequencies of 
rotation. 

54. The method of any one of claims 42 to 53 wherein the at least threei 
gradiometers are rotated at differing frequencies, in order to facilitate 
separation of their data signals in the Fourier domain. 

35 55. A device for obtaining a complete magnetic gradient tensor of a magnetic 
field, the device comprising: 
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at least three axial gradiometers positioned such that an axis of each 
axial gradiometer is not parallel to an axis of any other one of the at least three 
gradiometers; and 

means for controllably altering a position of each of the at least three 
5 gradiometers relative to the magnetic field during operation of the gradiometer. 

56. The device of claim 55 wherein the means for controllably altering are 
means for rotating each of the at least three axial gradiometers about a 
respective axis of rotation of that gradiometer. 

57. The device of claim 56 wherein the axis of rotation of each gradiometer 
10 is not parallel to an axis of rotation of any other one of the at least three 

gradiometers. 

58. The device of claim 56 or claim 57, wherein the at least throe 
gradiometers are rotated at a substantially constant velocity about their 
respective axes of rotation. 

15 59. The device of claim 58, further comprising means for converting data 
obtained from the at least three gradiometers to the Fourier domain. 
60. The method of claim 58 or claim 59 wherein the at least three 
gradiometers are rotated at differing frequencies, in order to facilitate 
separation of their data signals in the Fourier domain. 

20 61. The device of any one of claims 56 to 60 further comprising means for 
applying a least-squares best fit solution to obtained data. 

62. The device of any one of claims 56 to 61 further comprising means for 
detecting and measuring a DC offset. 

63. The device of claim 62 wherein the DC offset comprises one or both of: 
25 low frequency drift in at least one field vector sensor of the at least three 

gradiometers; and the fixed offset of at least one field vector sensor of the at. 
least three gradiometers. 

64. The device of claim 62 or claim 63 further comprising means for 
determining operating conditions of one or more of the at least three 

30 gradiometers based on said detected and measured DC offset. 

65. The device of any one of claims 62 to 64 wherein said means for 
detecting and measuring a DC offset separates the DC offset from field data 
and gradient data in the Fourier domain. 

66. The device of any one of claims 56 to 65 further comprising means for 
35 distinguishing field gradient information from field information in the Fourier 

domain. 
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67. The device of any one of claims 56 to 66 further comprising means for 
distinguishing information about the gxy component of the gradient tensor from 
information due to the diagonal gradient components, even at the same 
frequency. 

5 68. The device of claim 59 wherein the Fourier transform is implemented by 
means of a fast Fourier transform algorithm. 

69. The device of claim 59 wherein the Fourier transform is implemented by 
sampling discrete data points throughout each cycle of rotation. 

70. The device of claim 69 wherein eight discrete data points are sampled. 
io 71. The device of claim 69 or 70 wherein the discrete data points are 

equidistant throughout the cycle of rotation. 

72. The device of any one of claims 56 to 71 further comprising means for 
isolating the at least three gradiometers from physical vibration occurring in the 
vicinity of their respective fundamental and first harmonic frequencies of 
15 rotation. 
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